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Abstract. This paper considers a family of active scalar equations with transport velocities which are more 
^sj singular by a derivative of order /3 than the active scalar. We prove that the equations with < /3 < 2 are 

,-H Lipschitz ill-posed for regular initial data. On the contrary, when < /3 < 1 we show local well-posedness for 

f^ patch-type weak solutions. 
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1. INTRODUCTION 



In this paper we study a singularly modified version of the incompressible porous media equation. We 

investigate the implications for the local well-posedness of the equations by modifying, with a fractional 

derivative, the constitutive relation between the scalar density and the convecting divergence free velocity 

vector. Our analysis is motivated by recent work ||2|| where it is shown that for the surface quasi-geostrophic 

equation such a singular modification of the constitutive law for the velocity, quite surprisingly still yields a 

T^ locally well-posed problem. In contrast, for the singular active scalar equation discussed in this paper, local 

C^ well-posedness does not hold for smooth solutions, but it does hold for certain weak solutions. 

H The incompressible porous media (IPM) equation itself is derived via Darcy's law for the evolution of a 

'— ' flow in a porous medium. In two dimensions it is the active scalar equation for the density field p{x, t) 

<^ dtp + v-vp = ^, (1.1) 

O V--u = 0, (1.2) 

(^ where the incompressible velocity field v is computed from p and the pressure P via Darcy's law lITlfTSl 

O n 

^ ^^v = -VP-g{Q,p). (1.3) 

o 

pg Here p is the viscosity of the fluid, k is the permeability of the medium, and g is gravity. For the sake of 

T-H simplicity, we let /i = k = <? = 1. The equations are set in either R^ x [0, oo) or T^ x [0, oo). We observe 

j>! that V is determined from p by a singular integral operator: using incompressibility ( |1.2[ ), the pressure is 

obtained from the density as P = {—/S)^^dx2P, which combined with ( |1.3[ ) yields 



-A)-l9,,p-(o,p) = i^^i^l/5, (1.4) 

where R = {Ri , i?2) is the vector of Riesz transforms. There is a considerable body of literature concerning 
the IPM equation (i.e., ( |1.1[ ) and (L4l), which as well as describing an important physical process, gives a 



simple model that captures the non-local structure of variable density incompressible fluids lITllSlfTBll. 

Active scalar equations that arise in fluid dynamics present a challenging set of problems in PDE. Maybe 
the best example is the surface quasi-geostrophic equation (SQG), introduced in the mathematical literature 
in in . Here the operator relating the velocity and the active scalar is also of order zero (as is the case for 
IPM). The SQG equation reads 

dte + u-ve = {), (1.5) 

u = R^d. (1.6) 
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For regular initial data, similar results have been proved for IPM and SQG Emm, while for weak solutions 
one can find different outcomes ||6l \T5\ fTF], and for patch-type weak solutions the two systems present 
completely different behaviors [7l|TT|. The global existence of smooth solutions remains open for both the 
IPM and SQG equation. 

There is a significant difference between the SQG and IPM equations which we explore in this paper: 



the operator in ( |1.4| ) is even, while the analogous operator in ( |1.6| ) is odd. In a recent paper [2|, the authors 
investigate what happens in a modified version of the SQG equations, when a (fractional) derivative loss 
in the map relating the scalar field and the velocity is included, i.e., instead of (|1.6|) one has u = A^i?"* 



Here /3 > and A = (— A)^/^ is the Zygmund operator. It is shown in [2| that since the divergence-free 
velocity u is obtained from 6 by the Fourier multiplier with symbol ik^\k\^^^, which is odd with respect 
to k, one obtains a crucial commutator term in the energy estimates. It then follows that the equations are 
locally well-posed in Sobolev spaces H^ with s > 4. In this paper we consider a modified version of the 
IPM equation where the fractional derivative A^ is inserted in the constitutive law ( |1.4[ ). More precisely we 
study the singular incompressible porous media (SIPM) equation, which is given by 

dtp + v-Vp = 0, (1.7) 

V = - V(-A)-ia^2 A^p - (0, A^p) = R^RiK^p, (1.8) 



where < /3 < 2. We observe several significant features of the operator in ( |1.8| ). It is a pseudo-differential 
operator of order /3, which is inhomogenous with respect to the coordinates xi and X2- Furthermore, it is 
an even operator, in the sense that its Fourier multiplier symbol, given explicitly as —kik-^\k\l^~'^ , is even 
with respect to the vector k. These features of the SIPM constitutive law ( |1.8| ) lead to results that are in 
sharp contrast to those obtained for the singular SQG in |2]. Considering ( |1.7| )-( [L8] ) with < /3 < 2 
on the spatially periodic domain, we prove that for smooth initial data the equations are locally Lipschitz 
ill-posed in Sobolev spaces H^ with s > 2. In contrast, when < /3 < 1, we prove local well-posedness 
for some patch-type weak solutions of ( |1.7| )-( [L8] ), on the full space. This dichotomy is a reflection of the 
subtle structure of the constitutive law ( |1.8| ). The even nature of the symbol relating the active scalar and the 
drift velocity also proved crucial in showing that for such active scalar equations L°° weak solutions are not 
unique IHEl. 

In Section 2 we prove Lipschitz ill-posedness in Sobolev spaces for the SIPM equation ( |1.7| )-( [L8] ) with 
< /3 < 2, on T^ X [0, oo). The proof follows the lines of a similar result for another active scalar equation 
where the constitutive law is given by an even unbounded Fourier multiplier, namely the magnetogeostrophic 
equation (MG) studied in llT0l[T2l . We use the techniques of continued fractions to construct a sequence 
of eigenfunctions for the operator obtained from linearizing the SIPM equation about a particular steady 
state. These C°° smooth eigenfunctions have real unstable eigenvalues with arbitrarily large magnitudes. 
Once such eigenvalues are exhibited for the linearized equation the Lipschitz ill-posedness (in the sense that 
there is no solution semigroup that has Lipschitz dependence on the initial data) of the full nonlinear SIPM 
equation is proved using classical arguments (see, for example, [ 10, JL4J). To emphasize that the crucial 



features of the operator in ( |1.8[ ) are that it is even and unbounded, we prove the Lipschitz ill-posedness of a 
general class of active scalar equations satisfying these properties. 

In Section 3 we consider weak solutions of ( |1.7| )-( [L8| ) when < /3 < 1, on M^ x [0, oo). We study 
solutions evolving irom patch-type initial data 

^(^ n) - / ^' '"^ Dl = {xeR^:x2> /o(xi)} 
where p^ > p^ are constants, and fo{xi) is a smooth function. Patch-type solutions evolving from such 



initial data are a priori more singular, and they have to be considered in a weak sense (see Definition 3.1 1. 
In particular, the velocity at points (xi, /o(3;i)) diverges to infinity. However, the evolution of these weak 
solutions can be reduced to a contour dynamics equation for the free boundary f{xi,t), with /(xi, 0) = 
/o(xi). This permits us to follow the construction given in [i8J for the classical IPM equation (/3 = 0). We 
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prove that for initial data of the form ( |1.9[ ) the SIPM equation with < /3 < 1 has a unique local in time 
patch-type solution, with a corresponding smooth interface / € L°°(0, T; H^{M.)), for s > 4 and T > 0. 
In the appendix we prove an abstract result concerning continued fractions, which is used in Section[2] 

2. ILL-POSEDNESS FOR ACTIVE SCALARS WITH EVEN UNBOUNDED CONSTITUTIVE LAWS 

In order to obtain the Lipschitz ill-posedness for the SIPM equation, we first study the equation linearized 
about a certain steady state. We prove that the associated linear operator has unbounded unstable spectrum. 



which is the main ingredient in the proof of Theorem 2.4 In the last subsection we prove that ill-posedness 



in Sobolev spaces holds for a general class of active scalar equations for which the velocity is obtained from 
the scalar via an even unbounded Fourier multiplier. 

2.1. Linear ill-posedness in L^ for the SIPM equation for < /3 < 2. Consider a density given by 

p{x,t) =P{X2) 

for a general function p : T — ^ M. From ( |1.8| ) we compute the steady state velocity as 

^ = (0,A^p)-(0,A^p) = (0,0) 

and therefore we get a steady state of ( |1.7| )-( [L8] ). 

Associated to the steady state p, one may define the linear operator L obtained by linearizing the nonlinear 
term in ( |1.7| ) about this steady state, namely 

Lp = -v-'Vp-v-'Vp = -V2d2P = -Rlk'^p d2p. (2.1) 

Using the method of continued fractions, see also |l9l [lOl, we shall prove next that the operator L has a 
sequence of eigenvalues with positive real part, which diverge to oo. This in turn implies that the linearized 
SIPM equation 

dtp = Lp (2.2) 



is ill-posed from H^ i— )• LfL^., for any s > 0. The singular features ( |2.2[ ) shall be used in Section 2.3 
to show that the full, nonlinear SIPM equations are ill-posed in Sobolev spaces, by using a classical per- 
turbation argument (see, for instance [.14.1 ). The following lemma is the key ingredient of the ill-posedness 
result. 

Lemma 2.1. Fix an integer a > 1, let s > 0, and chose the steady state p{x2) = sin(ax2) of ( |1.7| )-( [L8] ). 
For any integer k > 1, the linear operator L associated to p, has a H^ smooth eigenfunction pkixi, X2), 
with WpkWH" = 1. and corresponding eigenvalue A^ > which satisfies 

for some constant Ca > 1, which is independent of k. 



Proof of Lemma 2.1 Following the arguments [10], we prove the lemma by explicitly constructing an eigen- 



function p, with associated eigenvalue A, i.e. a solution of 

Lp = -acos(ax2)i??A^/9 = Xp, (2.4) 

where we make the ansatz that p is given explicitly by the Fourier series 

p{xi,X2) = sin(A;xi) N^ c„ sm{nax2) , (2.5) 



n>l 



where c„ 7^ for all n > 1. Note that given A and ci, one may solve (2.4 1-(2.5 I for all Cn, with n > 2, but 



only for suitable values of A do these Cn's converge sufficiently fast to as n — > cx). 
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Inserting the ansatz ( |2.5| ) into ( |2.4[ ), and matciiing the terms with same oscillation frequency, one obtains 
that recursion relation 



Aci + 



C2 
P2 



. , C„+l C„_l 
XCn H ^ 



0, when n = 1, 

0, for all n> 2, 



Pn+l Pn-1 



where for all n > 1 we have denoted 



Pn 






(2.6) 
(2.7) 

(2.8) 



Note that pn grow unboundedly as n — )■ oo, whenever /3 < 2, and they are monotonically increasing. To 
solve (2^i-(2^i it is standard to introduce rjn = {cnPn-i) / (cn-iPn) , which solves 

Api + i]2 = 0, when n = 1, 



Xpn + Vn+l + 



1 



??»i 



0, for all n > 2. 



(2.9) 
(2.10) 



Note that if A is known, the recursion (2.9l-( 2.10 i may be used to determine the values of {r]n}n>2, and 
also the sequence {c„} by setting ci = pi and 



■ PnVn ■■ -m 



(2.11) 



for all n > 2. The compatibility of ( |2.9| ) and ( |2.10| ) requires that A is given by a root of the characteristic 
equation 

1 



Api 



Ap2 



(2.12) 



Ap3- 



Ap4- 



To solve the characteristic equation ( |2.12[ ) we appeal to the following abstract result about continuous frac- 
tions, whose proof we give in Appendix [A| below. 



Theorem 2.2. Assume that the sequence of real numbers {pn}n>i satisfies 

<Pn< Pn+l 

for any n > 1, and the Pn 's are unbounded, that is 



lim pr, 



oo. 



Then there exists a real positive solution A* of ( 2.12 ), such that 

< A=^ < 



y/PlP2 



\JP\P2 - P\ 



(2.13) 

(2.14) 
(2.15) 



holds, and the sequence {cn}n>2 defined by (2.11 1 decays exponentially fast for large enough n and we 
have the estimate 



l™''Cn||^2(p^) < C i |no|'' + 



1^2 (N) 



(2.16) 



where uq is the largest integer such that pn^ < 4p2, cind C > is a constant. 



Since the sequence {pn} defined in (2^1 above is monotonically increasing and unbounded when < 

3f 

1 



/3 < 2, we may apply Theorem 2.2 and obtain the existence of a solution A^ to ( |2.12| ), which satisfies 

1 



VPi^ 



< A^ < 



VPiP2 - Pi 



(2.17) 



SINGULAR INCOMPRESSIBLE POROUS MEDIA EQUATION 



Note that in addition to the existence of Afc, Theorem 2.2 also guarantees that the coefficients {cn} con- 
structed via ( |2.11| ) decay exponentially fast after large enough n, so that the function constructed in ( |2.5| ) is 
smooth, and in particular lies in H^. Moreover, letting 

n>l 

we may divide the Cn shy Cs,k, and define 

pixi,X2) 



Pkixi,X2) 



C. 



s,k 



which is still a smooth eigenfunction of L with eigenvalue A^, and is normalized to have unit H^ norm. 
Under this normalization, in view of (|2.16|) we may also estimate the L^ norm of p^ 



II II -^11 II 1 ||ra"Cn||^2(M) _ 1 ^ 1 

WPkU^ - cJ'^nWmn) > cinf,+p2/pi) Cs,k ~ C7K+P2M) - CaA^ + k^y ^ ^ 

where Ca,s is a positive constant which depends only on a and s. Above we have used that no is the largest 
number such that p„Q < 4pi, which in view of ( |2.8[ ) may be computed explicitly (no ~ 2k /a), and P2/P1 is 
uniformly bounded in k. 



To conclude the proof of the lemma it is only left to verify that \23\ holds. Inserting the exact form of 



p\ and p2 from ( |2.8| ) into the estimate ( |2.17[ ), yields the existence of a positive constant Ca such that ( |2.3| ) 



holds, thereby concluding the proof of the lemma. D 



2.2. Linear ill-posedness for the SIPM equation when /3 = 2. When /3 = 2 the argument in Section 2. 1 
does not apply directly since the corresponding pn defined by ( |2.8| ) are not growing unboundedly. However 
the linear ill-posedness still holds, as the equation is even more singular. In fact the argument given here 
works for /3 < 4. Here we do not construct an explicit eigenvalue of the linearized operator, but instead give 
a lower bound for the solution at time t. 

Again we linearize ( |1.7| )-( [L8] ) around p = sin(ax2). The associated linearized operator is Lp = 
— acos(ax2)-Rf A^yO. To show the linear instability of dtp = Lp, define another linear operator 

Lp = —acos{ax2)RiA.p 

which corresponds to the operator of ( |2.1[) w ith /3 = 1. Let {Xk,Pk) be the corresponding sequence of 



eigen-pairs of L as constructed in Section 2.1 Then there exists Ca such that kC^ ^ < A^ < Cak'^. By the 



definition of p^ in ( |2.5| ), we have fj2 Pk = 0. We now define 

Pk = A~ p/j = sin(A;xi) > — sin(nax2) , 



n>l 

Then 

Lpk = -acos{ax2)RlApk = ^kPk = >^k^Pk- 
Therefore, the unique solution of 

dtp = Lp, p{0,x) = pk{x) , 
is p{t, x) = e*^'^^pk{x), which in turn shows that 

\\p{;t)\\L2>e'^'^\\pk\\L2, 

since pk has zero mean on T^. Hence ||p(-, i)||^2 > e*^* ||/9(-, 0)11^2, and since Xk can be made arbitrarily 
large by sending A; — )• 00, it follows that the linearized equations are ill-posed in the L^ norm (in the sense 
that there is no continuous semigroup at t = 0). 
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2.3. Nonlinear ill-posedness in H'^ for tiie SIPM equations. We recall of. llTOll the definition of Lipschitz 
local well-posedness: 

Definition 2.3. Let Y C X C L^ be Banach spaces. The initial value problem for the SIPM equation ( |1.7[ )- 
( 1.8 1 is called locally Lipschitz {X, Y) well-posed, if there exist continuous functions T : [0, oo)^ — )■ (0, oo) 



non-increasing (with respect to both variables), and K : [0,oo)^ — )■ (0, oo) non-decreasing, so that for 
every pair of initial data Pq,Pq € Y there exist unique solutions p^^> , p^"^' € L°° (0, T;X) of the initial 
value problem associated to ( |1.7| )-( [L8] ), that satisfy 

||^(i)(.,i)_p{2)(.,i)||^< ^11^(1) _p(2)||^ (2.19) 

for every t G [0,r]. Here T = r(||pi')||y, WpfWv) and K = K{\\p^^\\y, ||/of ||y). 

The Banach spaces X, Y considered here are X = H^ and Y = H^, with r > and s > r + 1. Indeed, 
if r < s < r + 1, the Lipschitz {H^, H^) well-posedness of first order equations should in general not even 
be expected, due to the derivative loss in the non-linearity. The main theorem of this section is: 

Theorem 2.4. The SIPM equations, with < /3 < 2, are locally Lipschitz, {H"^ , H^) ill-posed, for any r > 2 



and s > r + 1, in the sense of Definition 2.3 above. 

The main idea of the proof of Theorem 2.4 is to let p^ be the steady state p{x2) = sin(aa;2), for some 



(2) _ 

fixed o > 1, and Pq = p + epk, where k is chosen to depend on the Lipschitz constant K and the time of 



existence T, in such a way that K < 2 e:K.p{k^T). Letting e — )■ it will follow that the linear equation ( |2.2| ) 
should be Lipschitz {X, L^) well-posed with the same Lipschitz constant K, on [0, T), which gives rise to 
a contradiction due to the choice of k. In order to implement this program we need to show uniqueness of 
solutions to the linearized SIPM equations. 

Proposition 2.5. Let p € L°^(0, T; L^(T^)) be a solution of the initial value problem 

dtp = Lp, p(-,0) = (2.20) 

where as before the linear operator L is defined as Lp = —a cos{ax2)R\^^ p, with a G Z, and < /3 < 2. 
Then for any t G (0, T) we have p{-,t) = 0. 

Proof of Proposition \2. 5\ We write p in terms of its Fourier series as p{x, t) = J2kez'^ p{^^ *) exp(ifc • x), 
and for each /ji G Z define 

[p{k,,t)f:=Y,\p{kiM,t)? 

fcaez 

which is finite for each fci G Z and t G (0, T) by the assumption p G L°°(0, T; L^). Taking the Fourier 
transform of ( |2.20 ), and using a cos(ax2) = a(e*"^2 + e~*°^2)/2, we obtain that 

dt[p{ki,t)f < a\kif Y^ {\p{ki,k2 - a,t)\ + \p{ki,k2 - a,t)\) \p{ki,k2,t)\ < 2a\kif[p{ki,t)f 
fc2ez 

for /3 < 2. The proof of the proposition is concluded since [p{ki, 0)] = for each ki G Z, as p{-, 0) = 0. □ 

Having established the uniqueness of solutions to the linearized equation, we now give the proof of the 
nonlinear ill-posedness result. 

Proof of Theorem [Z4] Fix throughout this proof a > 1 and p{x2) = sin(ax2) a steady state of ( |1.7[ )-( [L8] ). 
Since r > 2, by assumption H^' is continuously embedded in H^ (for any /3 G (0, 2]), and the linear operator 

Lp = -RlA^p d2P 

maps X continuously into L^. In addition, the nonlinearity 

Np = R^RiK^p ■ Vp 
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may be bounded as 

\\Np\\l2 < C\\A^p\\l2\\Vp\\l^ < C\\pfHr (2.21) 

for some constant C > 0, since in two dimensions H^^^ C L°° for r > 2. Assume ad absurdum that the 
SIPM equations are locally Lipschitz {X, Y) well-posed in the sense of Definition 2.3 

We let Pq (x) = p{x2), so that p^^'{x, t) = p{x2) is the unique solution in H''' of ( |1.7| )-( [L8| ) with initial 
data Pq . Denote ||p||h'' for simplicity by C. Let t/^q G H'^ be a smooth function, to be chosen precisely 
later, such that HV'oIIh^ = 1- For each e G (0, C] we may define p^ ' {x) = Pq{x) = p{x2) + e'4)Q{x) e H^, 
and we denote the unique solution in H^ of ( |1.7[ )-( [L8] ) with initial data pg by p^ (instead of p^^'^^). By 
Definition 2.3 there exits a time T^ = T^{C, \\Pq\\h^) and a Lipschitz constant iCg = K^{C, \\Pq\\h^), such 



that we have 



sup I 

since ||^o||h'= = 1- Note that | 
assumptions on the functions T{ 
constant K > such that we have 



\p'{-,t)--p{-)\\Hr<K, 



\PI 



P\\h^ 



K,e 



PoWh 



< IIpIIh" + £ ^ 2C for all e G (0,C], and hence due to our 
and K {•,-), there exists a time of existence T > and a Lipschitz 



sup 

[o,f] 



\p'{-^t)-p{-)\\Hr<Ke 



(2.22) 



for any e G (0, C]. That is, T^ and K^ may be chosen independently on e. 

In view of the definition of pg, we have that iI)q = (pq — p)/e, and we may write the solution p^ as an 
0{e) perturbation of p, i.e. 



It follows from ( |2.22[ ) that {V'^Je is uniformly bounded in L°°(0, T; ff ) by K and V'^ is a solution of 

OiV" = Li;, + eA^V'^ V"(-, 0) = ^o- (2.23) 

By ( |Z2T] ) we infer that 



WN^ 



|L2 



<c\ 



nlr 



< CK^ 



(2.24) 



on [0, T], and hence since i?^^ C H^^ we infer from ( |2.23"| ) that {dt4>^}e is uniformly bounded in L°°(0,T; L^) 
by CK^ + -^||(92p||l°°- Therefore, by the classical Aubin-Lions compactness lemma we obtain that the 
weak-* limit if) G L°^(0, T; W^) is such that ip^ ^ t/; strongly in the L^ norm. But sending e — )• in ( |2.23| ), 



by using ( |2.24| ) we obtain that ip is the unique solution of the initial value problem 

dt^ = Lllj, -(/;(., 0) = V-o, 

and satisfies 

sup||V'(-,i)||L2<^- 

[0,f] 



(2.25) 



(2.26) 



Uniqueness follows from Proposition 2.5 above, since ( |2.25| ) is a linear problem. 



The proof of the theorem is now concluded b y car efully choosing the initial data ^o G H^ of p.25| ), in 
terms of T and K. More precisely, by Lemma 2.1 for any A; > 1 we may find a smooth eigenfunction 
Pk of the operator L, normalized to have H'^ norm equal to 1, such that its associated eigenvalue satisfies 
Afc > k^ /Ca (where Ca is a positive constant). It follows that the solution '>p{x,t) of p.25| ) with initial 
condition -00 = Pfc> is given by e^p{tXk)ipo{x) (again we invoke Proposition 2.5 for uniqueness). Therefore, 
recalling how p^ was constructed, by ( |2.18| ) we obtain 



,r)||i2 = eXp(rAA;)||pfe||L2 > 



exp(f Afc) ^ exp{fk^/Ca) 



^a,sk 



> 



^a,sk 



(2.27) 
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where Ca and Ca,s are constant that may depend on a and s. Since for any given T,K > 0, we can find a 
sufficiently large k such that ex.p{Tk^ / Ca) / {Ca,sk^) > 2K the proof is now completed, since we arrived 
at a contraction with (12. 261). D 



2.4. lU-posedness of active scalar equations with singular even constitutive law. The method used to 
prove ill-posedness for the SIPM equations may be directly generalized to show the ill-posedness for a class 
of active scalar equations of the type 

dte + u-V9 = 0, (2.28) 

V • -a = 0, u = MO, (2.29) 

where {x, t) G T'^ x [0, oo). The d-dimensional vector field u is obtained from via the Fourier multiplier 
operator M, which is given explicitly in term of the Fourier symbol m = {mi, . . . , md-i,md) : Z"^ — >• M''. 
We denote the frequency variable by k. In this section we give sufficient conditions for M which ensure 
the ill-posedness of ( |2.28| )-( |2^ . 

Let j G {1, . . . , d} be a fixed coordinate, which for ease of notation we simply take to be j = d. We 
write k' to denote the d — 1 dimensional vector (/ci, . . . , kd-i) G If^^^. We assume the following hold: 

(i) m,(0',a) = for a given positive integer a (which we fixed throughout this section); that is, 
9 = sm.{axd) is a steady state solution of ( 2.28 l-( 2.29 ) with corresponding velocity u = Q; 

(ii) md{k) is a real positive rational function, that is even in k; 

(iii) md{k' , na) — )• cxd as \k'\ — ;• cxd, for any fixed n G N; 

(iv) md{k' , na) — )• as n — )• oo, for any fixed k' G Ij'^~^; 

(v) md{k' , {n + l)a) < md{k' , na) for all n G N, and any fixed k' G Z^^^; 

(vi) |m(fc)| < C(l + |/c|)'^o for some ro > and C > 0, for all k G Z*^. 

Examples of such equations are given by the magneto-geostrophic equation introduced in lfT2]| (see also 
11 Oil ), and the singular incompressible porous media equation, both in two and three dimensions. 

Theorem 2.6. Assume the Fourier multiplier symbol m satisfies properties ^—^ above. Then the active 
scalar equation ( 2.28 l-( 2.29 ) is Lipschitz {H''' , H^) ill-posed, for any r > max{2, tq} and s > r + 1, in 
the sense of Definition\2.3\ 



Proof. We will only prove the unboundedness of the spectrum of the linearized operator associated with 
( |2.28| )-( |2.29[ ). The nonlinear Lipschitz ill-posedness follows by arguments verbatim to those in Section [23] 
It is here where condition ([vi]l is used: to ensure that for r > max{ro, 2} the bound analogous to ( |2.21| ) 
holds. We omit these details. 

By assumption ^, 6{xd) = sin(axd) is a steady state solution of ( 2.28 1, for some given a G N. The 
linearized operator around 0{xd) has the form 



Le{x) = -Me{x) ■ Ve{xd) = -Md9{x) e'{xd) = -{mde)^{x)acos{axd) 
We will construct an appropriate eigen-pair (A, 9) of L, i.e., a solution of 

Le = \e. 

Fix k' G U^^'^, and make the ansatz 

d-l 
^(*) = II sm{kiXi) ^ Cn sm{naxd). 



(2.30) 
(2.31) 



(2.32) 



i=l 



n>l 



Inserting ( |2.32| ) into ( |2.31| ) and matching the corresponding Fourier modes gives the recursion relation 

C2 



Aci + 



P2 

XCn H \ 

Pn+1 Pn-1 



0, when n = 1, 

0, for all n > 2, 



(2.33) 
(2.34) 
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where for all n > 1 we have denoted 

Pn = f/° , . (2.35) 

mci[k ,naj 

We point out that by ([jv]l and ([v]) we know that the p„ are monotone increasing, and growing unboundedly. 
Hence, to solve ( |2.33[ )- p.34| ), as in Section 2. 1 we need to find a positive root of the characteristic equation 



Xpi = —^ . (2.36) 



This root exists in view of Theorem 2.2 since the p„'s are increasing and unbounded. We also obtain that the 



coefficients c„ decay exponentially fast so that the function 6 constructed in ( |2.32[ ) is smooth. In addition, 
we have a bound for A of the form 

' <A< ' 



Vp^ \/piP2 - pI 

which combined with dm]) shows that we can find arbitrarily large A, by simply letting |A;'| be large enough. 



This shows that the spectrum of the linearized operator L contains eigenvalues of arbitrary large positive 
part, and therefore the linearized equation is ill-posed, in the sense that it possesses no semigroup that is 
continuous at t = in L^. D 

3. Local well-posedness for weak solutions of patch-type for the SIPM equations 
In this section we consider solutions for a scalar p{x, t) given by 

,(^ A-i P' in D\t) = {x&^^:x2> f{xi,t)} 

where p^, p^ > are constants, p^ ^ p^ and the common boundary dD^{t) j = 1, 2 is parameterized as 
^2 = f{xi, t). If p{x, t) satisfies ( |3.1| ), ( |1.7[ ) and ( L8 1 we say that it is a patch-type solution. Then SIPM is 



understood in the distributional sense and its precise definition is as follows: 

Definition 3.1. Let T > 0. A function p G L°°{0, T; L°°(R'^)) satisfying ( |3TT] ) is a weak solution of ^TJ) - 
( |1.8| l if for any test function cp £ C^([0, T) x M?), the following integral equation holds: 

/ / p{dt4> + v-'V(f))dxdt+ po{x)(t>{x,0)dx = 0, (3.2) 

Jo iR2 iR2 

and V may be computed from p by means of (|1.8|), i.e. v = R RiA^p in the sense of distributions. 



We will show later that if p satisfies ( |3.1| ) and / is smooth then 



\v{x,t)\ < Cj^^^^j^^ e LL(M2),Vt > 0, (3.3) 

due to < /? < 1. This implies that the nonlinear term in ( |3.2[ ) is well defined and may be bounded as 

rT 



/ / \pv-V(l)\dxdt<\\p\\L^\\v\\Li(^^^pp^ 

Jo JR2 



Here we give the main ingredients to obtain the following contour equation for /: 

p'-p' f {ri-C){d,f{v,t)-d,f{C,t)) 



Cp A((7?-C)2 + (/(r/,t)-/(C,t))2)^ 
fir],0) = foir]), (3.5) 

where ry G M, C/3 > 0, and < /3 < 1. Next we obtain local-existence for the system above with p'^ > p^. 



10 



SUSAN FRIEDLANDER, FRANCISCO GANCEDO, WEIRAN SUN, AND VLAD VICOL 



To get the evolution for / we need the velocity at (r/, f{r], t)) but only in the normal direction. In fact 
iv, fiv, t))t ■ i-d^fiv, t),l)= v{r], f{rj, t), t) ■ {-d^fiv, t), 1). 
We understand above expression for the velocity with the following limit 



v{rj,f{r],t),t)-{-d^f{rj,t),l) = lhnv{r]-edrjf{v,t)J{v,t) + e)-{-dr,f{v-edrjf{v,t),t),l). 



(3.6) 



We shall now prove that the limit in ( |3.6| ) is exactly the expression on the right side of ( |3.4[ ). By ( |1.8[ ) we 
have the following relation between v and p 

V = a^,A-2+/3v^p. (3.7) 

For g a regular function, it is a classical fact that 

1 f xi-yi 



d^.A-'+fgix) 



Cb 



\x-y 



2+13 



g{y)dy 



(3.8) 



where Cp = (7r22~^r(^))/(/3r(/3/2)) is a normaUzation constant. The identity 

V^p{x, t) = (p2 _ pi)(i, d^f{7], t))6{x2 - f{ri, t)), 
where 6 stands of the Dirac delta function, combined with (|3.7|l and (13. Sb allows us to write 



-(x,t) = -^/i^i:i#%MCJ^c. 



Cb 



\x-{CJ{C,tW+^ 



(3.9) 



Using ( |3.9| ) we compute the limit in ( |3.6| ) 

,firj,t),fiv,t)^ 

P^-p' f {v-C-edr,firj,t)){dr,f{C,t)-dr,f{v-edrjf{ri,t),t))dC 



limv{r] - ednfiv, t), f{r], t) + e) ■ {-dnf{r] - ed^f{r], t),t), 1) 



lim - 

e-S-0 



C, 



/3 



iiv-C- edr,f{v, t)r + [fiv, t) - /(C, t) + e)2) 



2+/3 



We split the integrand L*^ into L|(t/, C) = -Zy'^(C)^)x(IC~^l > '") ^nd L|(r/, C) = -^"^(C) ^)x(IC~^l ^ ''')' 
where x stands for the characteristic function, and r > is a fixed number. Without loss of generality we 
may take e < r/(2||9^/||ioo), since we send it to anyway. For such e we bound L\ pointwise as 

ire. .^| ^ \r] - C - edrjf jr], t)\ \drjf{C,t)-d^firj-€dr,f{v,t),t)\ 

jX[\C - m > r) < -: -:TTiiTX[\C-m>r). 



Since the right side of ( 3.10 1 lies in L^ (M) , and is independent of e, from the dominated convergence theorem 
we obtain that 



lim / Lii,,c)dC = [ (^-C)(V(C,t)-V(,,t)K 

^-0 Jr J\n~C\>r ((7? - C)2 + (/(r/, t) - fiC, t))2) -^ 



(3.11) 



On the other hand, we now show that the integral of LI is small uniformly in e. By the mean value theorem, 
and the fact that e < r/(2||9^/||Loo), we have 



\L'2iV,C)\dC<C\\dr,^f\\i 
< C\\dr,^f\\l 



|C-r,|<r \v-C-^dnf{v,t)\l^ 



\z\<2r 



dz 



< C\\dr,r,f\\L^r 



1-/3 



(3.12) 
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By combining ( |3.11| ) and ( |3.12[ ) we conclude that for any r > 

{v-C){dr,fiC,t)-d,fir],t))dC 



lini 

e->0 



L^{v,C)dC 



2+13 



which converges to as r — )• 0, since /3 G (0, 1). We have thus proven that 

lim v{rj - ed^f{r], t),f{r], t) + e) ■ (-5,,/(r/ - ed^f{r], t),t), 1) 



< CWd^r^fUoor^-^, (3.13) 



P - P 



(^-C)(5,/(ry,t)-a,/(C,tM 



{{r}-CY + U{v,t)-f[C,t)Y)-^ 



2+/3 ' 



and hence ( |3.4[ ) holds. The rest of the section is devoted to proving the following result. 

Theorem 3.2. Let p^ > p^, (3 £ (0, 1), and /o G H^for s > 4. Then there exists T = T{\\fo\\H^) > such 
that the contour differential equation given by ( |3.4[ )-( [331 ) has a unique solution f G C([0, T], H'^). 



We give the proof for s = 4 and leave s > 4 to the reader. For notational convenience, we take the 
coefficient (p^ — p^)/Cis = 1 and omit the time dependence of /. 

Proof of Theorem \T2\ We proceed by proving an a priori energy estimate of the form 



d 
dt 



H^<C{l + \\f\\H^] 



for C and k > 1 universal constants. By ( |3.4[ ) we have 

1 d 

2dt 
where 



2 
L2 



fir])ft{v)dv = Ii + h, 



and 



f 'I 
/ 'I 

J\c\<i Jr 



{\c\' + {m-fiv-cw)^ 

f{v){dr,f{r])-d,f{7j-C)) 



2+13 



dr]d(. 



'ICI<1 J^ i\C\^ + {f{r])-f {7] -CW)^ 
Using the Cauchy-Schwartz inequality and /3 > 0, we estimate 

^ ^ \f{vmdr,fm + \dr,f{v - 0\)dvdC < CWfWl, 



h < 



ICI>i 



ICI 



1+/3 



On the other hand, using the Cauchy-Schwartz inequality and the Gagliardo characterization of the Sobolev 
norm ifTTl we obtain 



l2< 



< 



ICI<i 



L2 



\f{r))\\d,f{v)-d,f {71-01 



drjdC, 



ICI<i 



ICI 



1+/3 



1/2 



<C||/|| 



L2 



ICI<i 



|C|i+/5 

drj{7l)-d^f[il-0?dr^j dC 

\ 1/2 



\d^f{ri)-d^f{il-C)Y 



ICI 



2+2,3 



dridC] <C7||/||i2||/||^i+,. 



Consequently 



Id^ 

2dt 



L2 



<C||/|| 



H2- 



(3.14) 
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Next we estimate the H^ norm of /. Using the Leibniz rule we obtain 

1 d 
2dt' 
where 



\d^f\\Ut) = [ d^fir])d^ftirj)dv = Jo + Ji + J2 + J3, 

JR 



Jr 



d^^f{v)d, / ad',-'f{v)-d:^-'f{v-owJC + {m-f{v-c)) 



2+/3 



dCdr] (3.15) 



for < i < 3. We deal first with the most singular term Jq, which contains fifth order derivatives. Integrat- 
ing by parts we rewrite 



•^0 = / d f{r])dr, 



d'm 



ie + im-fiv-cm^ 

2+0 "S."^- 

(C2 + (/(^)-/(r?-C))2)^ 



By the change of variable (r/, C) 



(ry — (, —(), we can further rewrite Jq as 



d'fiv-C) 



2T0" "C"^ • 

{e + {f{ri)-f{v-CW)^ 



Taking the average of ( 3.16 l and ( 3.17 ), and applying a further change of variables (??, C) 
have 



(3.16) 



(3.17) 



(??+|,C),we 



Jn 



{C + {f{v)-firi-C)r) 2 



2+/3 



(C2 + (/(r/+i)-/(r/-|))2)^ 



9cKV(^?+i)-5'/(r?-i))- 



W(^ + |)-5,V(r? 



|))C 



2+/3 



dQdrj 



dQdrj. 



Here we used that 9^/(r/ + C/2) = 25^5^/(r? + C/2) and d^fii] - C/2) = -29^9^/(7/ - C/2). Integration 
by parts in C, and the change of variables (r/, C) ~^ (^ ~ C/2, C) then gives 

Jo = l[ [ (5,V(r? + i) - 5,V(r? - i))'5c (C(C' + ifiv + i) - /(^ - l))')-"^) dCdv 

= \j I (5,V(r?) - 9,V(r? - C))'5c (C(C' + (/(r/) - /(r? - 0?)-"^) dCdrj. 

We compute explicitly the derivative term and obtain 

dc{ae + {firi)-f{v-0?f-^) 
_ (2 + (3){f{rj) - f{rj - C))(/(r/) - f{v - C) " ^./(t? - C)C) 



(3.18) 



1 + /3 



(C2 + (/(r/)-/(r/-C))2)^ 
Applying the explicit form of the derivative in ( |3.18| ) gives 

Jo 
where 

,(i)_2 + /3 



ic + ifiv)- fir] -oyf-^' 



/^^ + j(') 



j: 



(8,v(,) - «:/(, - 0)=^^- ■"'' - ^^»'^'''* - f^" - ^' - ^•'^<'' - '^'^' 



4+;3 



(C2 + (/(r?)-/(r?-C))2)^ 



(3.19) 

dQdrj 
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and 



J, 



(2) 



1 + /3 



{e + {f{v)-f{^-cwf-^ 

(1) 



Using Gagliardo's characterization of the Sobolev norm, we may bound Jk as 



A'^ < ^ 



iotm - otfiv - c))^ i^(-^)-^(-^:^)ii,^(-^)-/(-^-^,)-i;/(-^-^)^i .c^. 



{e + if{v)-fiv-CW)^ 



<c^II/IIh/^-II/II 



H/2 



idtf{v)-dU-iv-C))' 



\c\ 



1+/3 



dQdrj 



C\\f\\w^,^\\f\\w^,^\\Aid'^f\\l,. 



Applying the Sobolev embedding H C W '°^, and interpolation, we obtain from (3.20i that 



1+/3 



7(1) 



4+2/3, 



j;,"<c\\f\\iv'ii + 



r(2) 



2+/3n , 
Hi ) + 



i|A^a,Vll i. 

8(l + ||5,/||i. 



2+;3 



(3.20) 



(3.21) 



Next we show that the first term of Jq in ( |3.19| ) gives a dissipation, which is in fact needed to close the 
energy estimate. To this end, we separate this term as follows: 

l + P [ [ {dtf{v)-dtf{r,-0? 



J, 



(2) 



2 



{e + {f{v)-f{v-cw)^ 

id',m-d',f{r,-C)r 



d(drj 



ICI 



2+/3 



2+p 



1 + 13 



2(l + ||5,/||i^) 
1 + /3 

2(l + ||5,/||i^) 



2+0 
2 



(l + ((/(r,)-/(r?-C))/C)2)^ 
(J^tM^^lfirj-Ol (2) ^ ^(2) 



||A^/94.||2 , x(2) .(2) 



(3.22) 



where 



J, 



(2) 
01 



1 + /3 



(a,V(r?) - a;^/(r? - c))^ 

|C|2+/3 

X ((1 + iifiv) - fiv - 0)IC?r^''-^^" - (1 + |5,/P)-(2+^)/2) dC^r,, 



and 



J, 



(2) 



1 + /3 



02 



(9,V(r/) - 5,V(r? - C))' 



|C|2+/3 



:i + ia,/p 



2+;3 



2+;3 



(i + ll9^/llioo)-— dci^. 



7(2) 



l+l 



Note that Jq2 is well-defined for / € H 2 , and we clearly have 



j(^) < 

Jg2 ^ U . 



(3.23) 



Thus we only need a bound for Jq^ . By a similar argument as in the estimates for Jq , we estimate 



(1 + ((/(^) - firi - C))K?)-^'^^^" - (1 + |5,/P)-(2+'^)/2 < Cll/ll 



\W^'' 



■ICI 



and therefore 



r(2) 



1+/3 



1 3+2/3/ 



J,ii^ <C||/||H.2,..||A25,V|li2 <C||/||^Y"(1 + 



2+. ^ (i + /?)||A^a,V||i, 



(3.24) 
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Combining (3.22 1, (3.23), and (3.24i, we obtain 






1 + /3 



wl|A'^5,V|li.+C7||/||if^(l + ||/||i+/) 



4(i + IIVIIi^) 

Inserting the bounds ( 3.21[ ) and ( 3.25| ) into ( |3.19 l, we finally obtain that 



Jn< 



1 + /3 



:i + ii5,/iii^) 



wiiA'^5,viii. + ii/iiS''^(i + ii/iiS-/) 



(3.25) 



(3.26) 



Estimating Jj, with 1 < i < 3, is direct since in these terms the higher order partial derivative in / one 
can find is 4. We may bound 



1+/3 



Y.J^<C\\ff^,'^{l + 



3^^^ ^ {i+m^^ omh 



m 



2+/3 



i=l 



(i + IIVIIi^)~ 



which together with the estimate for Jq arising from (3.26l gives 



1^ 
2dt 



1+/3 



3+2,9, 



|9,Vlli2<C7||/||^r'^(l + 



3+. (i + /3)l|A^a^/lli. 

//4 J 



2+/3 



(1+115,/llioo)^ 



Finally, using the L estimate (3.14i we obtain 



2dt 



H^ 



<C\\f&{l + 



3+. (i + /3)l|A^a,Vlli, 

Hi ) 



1 



|9,/||i. 



2+/3 



(3.27) 



This concludes the a priori estimates needed to obtain the local-existence of smooth solutions to ( |3.4| )-( |33] ). 
Regarding uniqueness, consider two solutions /i and /2 of (34 i-(3^l. Using the nonlinear dissipation, 
a similar approach allows us to get 

ld„„ „„o „„ „ „„o (l+/3)||A^(/i-/2)||i, 



2dt 



l/l - /2||l2 < ^(11/111^4, 11/211^4)11/1 - /2|Il2 



2+,9 
I 2 



8(i + ||a^/illi. 

for some polynomial C(-, •). The Gronwall inequality then concludes the proof of uniqueness of solutions 
to (|X4l)-(|331l. 



The construction of solutions obeying the a priori estimate (3.27i, and verifying that the solutions to 



( |3.4[ )-( |331 ) give a weak solution of patch-type to the SIPM equations in the sense of Definition 3.1 follows 
precisely the same arguments as in the case of the classical porous media equation (i.e. /3 = 0). We omit 
these details here and refer the reader to [71|8l. D 



Appendix A. Proof of abstract results about continued fractions 



The goal of this appendix is to give a proof of Theorem 2.2 Recall that {pn\n>i is a strictly increas- 
ing sequence of positive numbers, and we need to find and estimate solutions A of the continued fraction 
equation 

1 



\pi 



Ap2 



Ap3- 



Ap4- 



which are real and positive. In addition, given such a solution A of ( |A. 1| ) we recursively solve 

ry„+i = — Apn for all n > 2 



r]n 



with initial condition 7^2 



for all n > 2, and ci = pi. 



-Xpr 

-Xpi, and then denote 

Cn = PnVn • • • % 



(A.l) 

(A.2) 
(A.3) 
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Proof of Theorem [Z2] For n > 2, and A > 2/p„ one may define the function 

G„(A) = '""-^fA^^ = y^ <A.4) 

Formally, we may write 

G„(A) 



A^r, 



\pn 



\pn — ■■■ 

for any A G V{Gn) = [2/pn,oo). Henceforth V will stand for domain. Note that V{Gn) C V{Gn+i) 
for any n > 2. Let observe some properties of the G^'s, that are all due in view of p.l3| ): G„(A) > 0; 
G„,(A) ^ as A ^ oo; l/(Ap„) < G„,(A) < 2/(Ap„); and G„+i(A) < G„(A) for any A G P(G„). Define 

1 2 

Ao = , Ai = — , and Vq = [minjAo, Ai}, oo). 

^/PlP2 P2 

One property of the Gn's which we will use later is that 

G„+i(A) < Xpn for any A G P(G„+i) n Vq (A.5) 

for any n > 2. Indeed, for those n such that Pn+i < '^Pn, we have 

2 
G„+i(A) < < Xpn, for any A G P(G„+i) , 

Xpn+l 

since p„,+i < 2pn =^ 2 < X^pnPn+i when A > 2/pn+i- If instead p„+i > 2p„, then for any A G Vq, 

1 4 K . r^ 



A PnPn+i > mill <^ PnPn+1, -^VnPn+i \ > min <^ — -, 4 ^ > 2 . 

\PIP2 pi J IPiP2 



which proves ( |A.5| ). 

For n > 2 we now introduce an auxiliary function F„(A) formally defined as 

i^n(A) = —, (A.6) 

■^Pn- -VZ 1 

^P"+l-Ap„+2-... 

for any A G Pq- In order to define F^ rigorously, for n > 2 and k > 0, define the (real) rational function 

Fn,k{X) 



Xpn 



>^Pn+l-- 



for any A > 2/pn-^-k+l■ After a certain value of k we have that all the rational functions F„ ^ are defined on 
Vq modulo a finite set of points where there are vertical asymptotes. Then we define 

F„(A) = lim F„,fc(A) 

for all A in Vq where this limit exits. At points where this limit does not exits, F„ has vertical asymptotes. 
Note that the equation we wish to solve, ( |A.1[ ) may now be written in this language as 

Xpi = F2iX). (A.7) 



In addition, if A* is a solution of (A.7 ) above, we have 



i^n+i(A*) = A*p„- — — — (A.8) 



for n > 2. Therefore Fn(A=„) = — r/„, where //„ was defined in ( |A.2| ) above. 
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We now show that Fn{X) is well-defined, continuous, and non-increasing for A G V{Gn) n Vq. For such 
a fixed A and n > 2, we have by ( |A.5 l that Xpn+k — Gn+k+i{X) > for any A; > 0. By ( |A.4| ) we also have 
Xpn+k - Gn+kW > and recalling Gn+k+iW < Gn+kW leads to 

< T ^ 7TV < T ^^ 7TV = G„+fc(A). (A.9) 



By iterating ( |A.9| ), in the and recalling the definition of F„ ^ one may show that 

< F„,fc+i(A) < F„,fc(A) < G„(A) 
holds for any A; > 0, and any A G V{Gn) n Vq. Due to the positivity of the i<^n,,fc's and using that 

Fn,k{X) 



Xpn — Fn+l,k-l{X) 

we in fact additionally obtain that l/(Ap„) < F„jc(A) for all k > 0. At last, we notice that since Gn+k+iW 
is decreasing in A, one may show that -F„ ^(A) is decreasing as well. Collecting the above obtained infor- 
mation of -Fri,fc. we may hence conclude that the function -F„(A) is well-defined, continuous, non-increasing 
and satisfies the bound 

-^ < F„(A) < Gn{X) < -^ (A. 10) 

Xpn Xpn 

ontheset2?(G„)nPo- 

In order to solve Api = F2(A) we distinguish two cases, depending on the relative size of Aq and Ai. The 
direct case is Ai < Aq, equivalently p2 > 4pi. In this case, inserting n = 2 in ( A.10| ) yields 

-^ < F2(A) < G2(A) 

Ap2 

for A > 2/p2 = Ai; thus F2 is continuous on "Dq- Due to the above inequality it is natural to check where 
Xpi intersects l/(Ap2) and G2(A). The graph of Xpi intersects the graph of l/(Ap2) at Aq = l/-^piP2 



and the graph of G2(A) at l/\/piP2 — pf- By the intermediate value theorem, there exists J—- < A* < 



Xj \Jp\P2 — p\ such that A*pi = F2(A*). When Ai < Aq we have thus found a solution A* to ( |A.1| ). 
For the A* found above, it is left to estimate the coefficients |c„| as defined in ( |A.3| ), i.e. 

|C„| =PnFn{X^)---F2(X*). 

But since A* > Ai = 1jp2, we may use the bound ( |A.10| ), and therefore 

2"-i pV.-^ 

|Cn| < —^. < — • 

A!^ 'Pn-l...p2 Pn-1.--P2 

Since the p„s are grow unboundedly, there exits an no > 2 such that pn^ < 2p2 < Pno+i- Due to the 
monotonicity of the p„s, we may therefore bound 

\Cn\ <P2 (A.ll) 

for all n < riQ and 

\Cn\ < T < — T < P22"°+^"'" (A.12) 

Pn-1 ■ ■ ■ Pno+lP2 Pno+1 



for all n > no + 1, which shows that the c„s eventually decay exponentially. Moreover, from ( |A.11[ ) and 
(lA. 12|l we obtain 



"0 / 2 \ 

n=0 fc>0 ^ ^1^ 



which proves ( |2.16| ). Here we used that ci= pi= P2{pi/p2)- 
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I G2(A) 




Figure 1 . Solving for A* when 2/p2 < l/\/PiP2- The above plot was obtained by numer- 



ically computing F2, where the coefficients p„ 



n 



arise in the study of the MG equation. 



In the numerical computation we have used the explicit formula for pn cf. |[TOl equation 
(2.33), and set all physical parameters to 1. 



The case Aq < Ai, or equivalently p2 < 4pi is more involved, since it will not be possible to show in 
general that F2 is continuous on Pq = [^^O; 00); we only know this on [Ai, 00). To overcome this we shall 



estimate the largest value of A where F2 has an asymptote, and work to the right of it. Due to ( |2.14[ ) there 
exits no > 2 such that 

2 2 

< Ao < — < Ai. 

Pno+l Pno 

Hence, by ( |A.10[ ) we have that F„g+i(A) is well-defined, continuous, and non-increasing on Vq = [Aq, 00). 



Moreover, by combining ( |A.10[ ) with ( |A.5| ) we see that 

< Fno+l{X) < Xpno 

on Vq. This hence allows to define the function F^g (A) on all of Vq, by setting 



\p. 



no 



no+1 



(A) 



(A. 13) 



Note that ( A.10[ ) only gives us that Fn^ is continuous, positive, and non-increasing on [2/p„(, , 00), but in 
view of ( A.13 l we now know these properties for Fno on all of Vq. Our goal is to iterate this process and 
define F„o_i, . . . , F2 on some large enough set. To achieve this we need to stay to the "right" of vertical 
asymptotes. We inductively define the sets 

2 



A^- = {A G Po: < Fj+i(A) < Xpj] D [- 
and on Aj we define the continuous, non-decreasing function 

F,(A) 



Pj 



,00 



APj - F'j+iW 



(A. 14) 



(A. 15) 



for all j G {2, no}. The induction starts at j = no, a case which was described in ( |A.13| ) and the paragraph 



below it. Intuitively speaking, if gj = inf ^j > Ao, since Fj+i is non-increasing, as A — ;■ aj+, we 
have that Fj{X) — > +00, i.e. Uj is the largest vertical asymptote of Fj. We observe that Aj is connected 
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(that is, an interval) since Fj+i is monotone and continuous on Aj. We either have that Aj 



Vq when 

Qj = Xq, or Aj = {aj, oo), with Oj E {Xo,2/pj). Note that the monotonicity of Fj follows from the 
monotonicity of Fj+i, which holds by induction, and from ( |A.15| ). We also note that by construction we 
have A2 C . . . C An^, = Vq; indeed if A G ^j-i, then < Fj{X) < Xpj^i, and hence by ( |A.15 l we must 
have < Fj_|_i(A) < Xpj, and so A G Aj. Hence, we finally obtain that 



no 



A, 



i=2 



is an interval which either equals Vq, or it equals (02, 00) for some 02 G (Aq, 2/^2)- Since -F3(A) > on 
A3D A2, by ( |A.15| ) we obtain that 

-^ < F2(A) (A. 16) 

Xp2 

on ^2- At last due to the monotonicity and continuity of F2 on A2, we obtain that the graph of Xpi must 
intersect the graph of F2(A) at some A=k G ^2- Indeed, if A2 = (02, 00), 02 > Aq, then limA-j>a2+ F2W = 
+00 and limA-j>oo ^2 (A) = 0, so we obtain the desired intersection point A=k > 02 from the intermediate 
value theorem. Otherwise, if A2 = Vq = [Ao,oo), we note that Xpi = l/{Xp2) at A = Ai, and so A^, 
such that A=kPi = F2(A*) must exist by ( |A.16| ) and the intermediate value theorem. In either case we have 
obtained a A=k > Aq that solves ( |A.1| ). In terms of upper bounds, either we have A* < Ai = 2/p2, or else we 
use that A* G [2/p2, 00) = T>{G2), and here we have F2(A) < G2(A) and hence A* < l/y/piP2 — pf, the 
intersection point of Xpi with G2 (A). Note that the later case can only occur if j?2 > 2pi, and in this case we 
can further estimate A* < l/pi- Thus, in general, we have obtained the upper bound A* < 1/ \/piP2 
stated in the theorem, and also the bound A,,, < max{2/p2, 1/pi}- 



Pi 




\/v\lv: 



Figure 2. Solving for A* when l/^/piP2 < 2/p2. and when F2 has vertical asymptotes 
larger than Aq. The above plot was obtained by numerically computing F2, where the 
coefficients p„ arise in the study of the SIPM equation. In the numerical computation we 
have used the explicit formula for p„ given in (|2.8[), and set all a = 6 = 1, and /3 = 1.5. 



In summary, we have proven the existence of a solution A* to ( |A.1| ), which satisfies the bound ( |2.15| ) 
Recalling that FniX^ = — ??n. as defined by ( |A.2| ), we obtain that the c„'s defined by ( |A.3 ) satisfy 



:p„i^„(A,)...F2(A,). 
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In order to estimate the | c„ | , we recall that one can find uq > 2 such that 

Pno < 4VPiP2 <Pno+l, 

or equivalently 2/p„„+i < Ao/2 < 2/pnQ. For this fixed no, using that X^, e Aj for all 2 < j < no, by 
( |A.14| ) we can estimate 

/ A\ \no-l 
Fno{X*)---F2{X,)<X::'-^Pno-l---P2< i^j (A.17) 

so that we have 

\cn\ < 4"n:!«-iAo"" < 4">2(A*v^I^)""-^ < P223"«-^ (A.18) 

for n < no, where we have also used the sharp bound A* < max{2/p2, 1/pi}, and p2 < 4pi. Then, using 
that for n > no + 1 we may use the bound ( |A.10| ), from the monotonicity of the pnS, and the choice of no 
we obtain from ( |A.17| ) that 

\cn\ = p„F„(A.)F„_i(A,) . . . F„„+i(A.)(4A.)"°-UJ-"° 

on— no 

< PnT^, (4A.)--iA^"o 

A* Pn ■ ■ ■ Pno+l 

^riQ — Ion— no ono / \ \ "~2no nSno— n 

< —tA —^ T < ^ 4^ < —. < P22^"°-" (A. 19) 

- ^n+l-2no^no-lpn-no-l - Ao V2AJ " Ao " ^' 

by recalling that by construction we have Aq < A*. We obtain that the c„'s decay exponentially for all 
n > 3no. Moreover, from ( A.18 1 and ( A.19 l we obtain 

3no 
Wn^r W'^n < r'^J7n^2s V^ 2 , 2 V^r,-2A: .^ ^ 2s|| ||2 _l9„2^^2s|| ||2 

IP Cn||^2(pj) ^ (,onoj 2^ '^n ^ P2 2_^ '^ ^ Ly?tO lFn||£2(j^) i" ^^2 ^ '-^'^0 l|Cn||^2(-p^), 

71=0 fc>0 



which proves ( |2. 16"1 ), and concludes the proof of the theorem. Here we also used p2 = C2/(piA*) < 2c2. D 
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